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ABSTRACT 

Many hot subdwarf stars show composite spectral energy distributions indicative of cool main 
sequence companions. Binary population synthesis (BPS) models demonstrate such systems can be 
formed via Roche lobe overflow or common envelope evolution but disagree on whether the resulting 
orbital periods will be long (years) or short (days). Few studies have been carried out to assess the 
orbital parameters of these spectroscopic composite binaries; current observations suggest the periods 
are long. To help address this problem, we selected fifteen moderately-bright (V~13) hot subdwarfs 
with F-K dwarf companions and monitored their radial velocities (RVs) from January 2005 to July 
2008 using the bench-mounted Medium Resolution Spectrograph on the Hobby-Eberly Telescope 
(HET). Here we describe the details of our observing, reduction, and analysis techniques and present 
preliminary results for all targets. By combining the HET data with recent observations from the 
Mercator telescope, we are able to calculate precise orbital solutions for three systems using more 
than 6 years of observations. We also present an up-to-date period histogram for all known hot 
subdwarf binaries, which suggests those with F-K main sequence companions tend to have orbital 
periods on the order of several years. Such long periods challenge the predictions of conventional 
BPS models, although a larger sample is needed for a thorough assessment of the models' predictive 
success. Lastly, one of our targets has an eccentric orbit, implying some composite-spectrum systems 
might have formerly been hierarchical triple systems, in which the inner binary merged to create the 
hot subdwarf. 

Subject headings: binaries: spectroscopic — ephemerides — subdwarfs — techniques: radial velocities 



1. INTRODUCTION 

Binarity plays an essential role in the story of the hot 
subdwarf B (sdB) stars, one of the most enigmatic stages 
of stellar evolution. The sdBs dominate surveys of faint 
blue objects and are found in almost all Galactic stellar 
populations. They are the field counterparts of the Ex- 
treme Horizontal Branch stars in globular clusters, and 
their location in the H-R diagram corresponds to stel- 
lar models with He-bur ning cores and extremely thin H 
envelopes (Heber 1986). Presumably, the sdB progeni- 
tors were stripped of almost all their surface hydrogen 
while on or near the red giant branch (RGB), igniting 
helium burning in their cores at nearly the same time. 
The question is how this can happen, since most Popula- 
tion I stars do not lose their entire envelopes, and remain 
in the red clump, not far from the RGB, during core He 
burning. 

Sto chastic mass loss from single stars (|D'Cruz et al.l 
1996) or other single-star scenarios to explain this are 
ad hoc, so most recent effort has concentrated on mass 
loss resulting from interactions in a clos e binary star 
system , following the ear l y exp loration by iMengel et al.1 
(|1976fl . lHan et al.1 (j2002t I2003T ) . for example, described 
five channels of binary star evolution that can lead to 
the formation of hot subdwarfs: the 'first' and 'second' 

* Based on observations obtained with the Hobby— Eberly Tele- 
scope, which is a joint project of the University of Texas at 
Austin, the Pennsylvania State University, Stanford University, 
Ludwig— Maximilians-Universitat Miinchen, and Georg-August- 
Universitat Gottingen. 
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stable Roche-lobe overflow (RLOF) scenarios, the 'first' 
and 'second' common envelope (CE) channels, and the 
merger of two He-core white dwarfs (WDs). Their CE 
scenarios produce close binaries with orbital periods from 
hours to days; the sdB companions are generally late- 
type (G/K/M) main sequence stars or white dwarfs. In 
contrast, their stable RLOF channels lead to systems 
with early- type (B/A/F) main sequence or late-type gi- 
ant companions and long orbital periods (P ~ few x 100 
d). An alternative sc enario employing the 7-form alism 
has been put for th by iNelemans et all (|2000l 120011 ) and 
INelemansT ([2010D . whose CE channels can produce sys- 
tems with main sequence companions and periods on the 
order of years. 

To date, orbital parameters have been measured for 
approximately 150 sdB binary systems; the overwhelm- 
ing majority of these are close binaries with P < 
15 d (see Table A.l of IGeier et all 120111 for a recent 
summary). The companions are mostly WDs and M 
dwarfs, essentially invisible in the glare of the sdB. 
With a suitable choice of tuning parameters, theoreti- 
cal models can reproduce the observed distribution of 
sdB mass es and orbital periods in these systems reason - 
ably well (|Maxted et allboOltlCopperwheat et alll2011l) . 
The recent binary p opulation synthesis (BPS) study by 
IClausen et all (|2012D . however, shows that numerous pa- 
rameter sets can reproduce this observed subpopulation; 
thus, additional constraints are needed to better tune BPS 
codes. 

Very little attention has been paid to the sdB+F/G/K 
binaries, in spite of the nearly equal flux contributions 
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TABLE 1 
RV Monitoring Targets 



Target 



Sp. Type" 


RA 


Doc 


V 




[J2000] 


[J2000] 




sdB+G2V 


00:42:06.1 


+05:09:24 


12.9 


sdB+GOV 


01:13:14.9 


+26:27:31 


12.9 


sdOV b +F4V 


01:23:43.2 


-05:05:45 


12.3 


sdB+G7V 


01:38:27.1 


+03:39:39 


13.4 


sdB+GlV 


02:35:12.0 


+09:45:38 


12.5 


sdB+GlV 


10:43:39.3 


+23:09:07 


13.4 


sdOB+G2V 


11:07:26.3 


+24:03:12 


11.3 


sdB+G9V c 


12:09:16.7 


+ 16:11:56 


13.8 


sdB+F8V 


12:56:04.9 


+28:07:20 


12.7 


sdO+GlV 


13:19:53.6 


+ 12:03:59 


11.3 


sdB+G4V 


13:40:14.7 


+60:52:48 


11.4 


sdB+G5 


14:50:36.1 


+65:05:53 


13.6 


sdOB+K5V c 


16:32:01.4 


+07:59:40 


12.8 


sdB+G8 


17:03:21.6 


+35:48:49 


13.2 


sdB+G4V 


17:19:45.5 


+51:52:10 


13.7 



Alternate Comments 
Name 



PG 0039+049 
PG 0110+262 
PB 8783 
PHL 1079 
PG 0232+095 
PG 1040+234 
PG 1104+243 
PG 1206+165 
PG 1253+284 
PG 1317+123 
PG 1338+611 
PG 1449+653 
PG 1629+081 
PG 1701+359 
PG 1718+519 



PB 6107 



EO Cct 



TON 1281 

PB 3854 
TON 139 
Feige 80 
Feige 87 



pulsating hot subdwarf 1 
resolved 1 

MS orbit published 4 



sdB orbit published' 
MS orbit published 4 
MS orbit published 4 

resolved 2 

resolved 1 



resolved ; presumed triple system 



"classifications from this work (sec i|4.H unless otherwise noted 
b sdO status of the hot subdwarf recently confirmed by^s^nsen (2^V£ 
"companion classification from Stark (2005) 



References: 1 Heber et al. (2002); 2 0stensen et al. (2005) : 3 Copperwheat et al. (2011) ; 4 0stensen U Van Winckel (2012) ; 5 Koen et al. 



from the two components in such systems. The relative 
dearth of published orbital parameters impli es their peri- 
ods ar e quite long. More than a decade ago. lGreen et al.l 
([20011 ) hinted that sdB binaries with composite spec- 
tra tend to have longer orbital periods. They reported 
an average around 3-4 years, although no specific sys- 
tems were named. Since that brief discussion, only two 
other studies have claimed the detection of orbital peri- 
ods in composite-spect rum binaries, both of which were 
published very recently. I0stensen fe Van Wincke] ()2012f ) 
(hereafter 012) discuss a survey conducted with the High 
Efficiency and Resolution Mercator Echelle Spectrograph 
(HERMES, see lRaskin et al.ll2011[ ) mounted at the Mer- 
cator Telescope on La Palma, Spain and present pre- 
liminary results showing orbital periods longer than ~ 
500 d for at least eigh t binaries. In addition to that 
work, iDeca et al.l (|2012l) report a period of 760 d for PG 
1018-047. 

Orbits with such long periods are sufficiently large to 
have once accommodated an inner binary, so it is possible 
that these systems were formerly hierarchical triple sys- 
tems, in which case the subdwar f could have form e d from 
the merger of the inner binary. IClausen fe Wadel (l20Tl . 
for example, recently proposed a scenario in which the 
merger of a He-core WD and a low-mass main sequence 
star could eventually form an sdB. If such a binary were 
in orbit with a dwarf companion, the merger would leave 
behind an sdB and a main sequence star that had no 
part in the formation of the sdB (aside from potentially 
advancing the merger process via the Kozai mechanism). 
Unlike stable RLOF- and CE-produced systems, which 
should have nearly circular orbits, no limitations exist on 
the eccentricities of these binaries other than a require- 
ment that the periastron separation not be too small. 
Thus, it is imperative to constrain the orbital geometry 
in addition to the period and radial velocity (RV) ampli- 
tudes. 

Here we describe a survey carried out with the Hobby- 
Eberly Telescope (HET) over w 3.5 years to measure the 
line-of-sight accelerations of a sample of sdBs with com- 
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posite spectra. We discuss our target selection in ij2]and 
observing strategies in £j3j The details of our analysis 
techniques for determining velocities and fitting the RV 
curves are given in Preliminary results and the cur- 
rent status of our program are presented in fj5j After 
combining the HET data with recent measurements from 
012, we report precise orbital solutions for three of the 
targets and discuss them in greater detail. Our initial re- 
sults confirm early suggestions that the orbital periods of 
many sdB+F/G/K systems are long, on the order of sev- 
eral years. We present in S}6] an updated orbital period 
histogram for all measured sdB binaries, which shows 
a possible gap in the period distribution centered near 
P = 70 days, and we briefly discuss the puzzle that it 
presents and possible resolutions. Finally, we summarize 
our results in fJTj 

2. TARGET SELECTION 

Our criteria for selecting targets were simple: tar- 
gets were required to be (i) hot subdwarfs with late F, 
G, or early K dwarf companions, (ii) relatively bright 
(V« 11-14), and (iii) observable from the HET, situated 
at McDonald Observatory. F/G/K main sequence com- 
panions have similar enough luminosities to hot subd- 
warfs that these systems show composite spectra at op- 
tical and/or infrared wavelengths. The cool compan- 
ions' rich spectra provide a plethora of lines for pre- 
cise RV measurement and match the solar spectrum 
(G2V) well enough that twilight sky spectra obtained 
nightly can be used as high signal-to-noise (S/N) ratio 
RV templates. Inf ormed by previous observations from 
IGreen et all (|200lb . we selected fifteen binaries meeting 
the above criteria, as presented in Table Q] All targets 
were recognized as poten tial hot subdwarf sy stems in the 
Palomar- Green Survey ()Green et alJ 1 1986^ . ex cept for 
PHL 1079 (lHaro fe Luvtedll962t lKilkennvlll984fl and PB 
8783 (jBerger fe Fringantl Il984f) . Rough spectral types 
for the cool companions have been published for most of 
the targets, but since these come from different sources 
and sometimes differ significantly, we classified them our- 
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selves as described in 34.11 

Several of th e targets have been obse rved in greater 
detail already. iHeber et al.l (|2002D and I0stensen et al.l 
(2005), for instance, found that four of the systems are 
resolved binaries with separations up to 2". Of these, 
PG 1253+284 has been recognized as a potential triple- 
star system due to the detection of a rapid acceleration of 
the sdB star that is incompatible with the separation im- 
plied by the 0".32 angular distance between the two ob- 
served component s. PB 8783 exhibits stellar pulsations 
(jKoen et al.lll997l ) with surface velocities that complicate 
the characterization of any orbital accelerations around 
its F-type companion. Most notably, three of our targets 
(PG 1104+243, PG 1317+123, and PG 1338+061) were 
recently observed in an RV monitoring program similar 
to ours by 012, who report circular orbit solutions from 
measurements of the cool companion lines. In fj5l we 
combine the HET data with their more recent results to 
gain a total observing baseline of ~6 yr for these binary 
systems. 

3. HET OBSERVATIONS AND REDUCTIONS 

Over a three-year period from 2005 to 2008, we 
monitored all survey targets spectroscopically with 
the be nch-mounted Medi um Resolution Spectrograph 
(MRS: lRamsev et al.lfT998h on the 9.2-m HET. The HET 
operates in a queue-scheduled mode. Follow-up observa- 
tions of some of our targets are currently being carried 
out with the High Resolution Spectrograph (HRS) on 
HET and will be discussed in a future paper. We used 
the "blue" 1.5" optical fiber pair, the 220 mm -1 cross- 
disperser grating, and 2x2 on-chip binning to achieve 
an average resolution of i?=10,000 over the 4400-6200 
A spectral range. This configuration yields an average 
dispersion of ~ 10 km s _1 per binned pixel and samples 
~3.2 pixels per resolution element. The pair of fibers, 
which are separated by 10", allowed us to monitor the 
sky flux and target flux simultaneously. Exposure times 
ranged from 180 to 1800 s depending on stellar magni- 
tude, giving us a typical S/N ratio of 80 per resolution 
element. 

Each object spectrum was preceded and/or followed 
by a pair of ThAr comparison lamp spectra (with dif- 
ferent exposure times) for wavelength calibration. Stan- 
dard calibration frames were collected each night, includ- 
ing twilight sky spectra, bias frames, and quartz lamp 
flat -field spectra. On most nights we also observed tel- 
luric, RV, and spectrophotometric standard stars using 
the same instrumental setup. 

We used the ccdproc routine in IRAE0 to bias-subtract 
and flat -field the spectra, then optimally extracted each 
aperture with the apall function. Finally, we used the 
ecidentify routine in iraf to identify comparison lamp 
emission lines and construct wavelength solutions. The 
ThAr solutions typically used ~ 100 lines and resulted in 
an RMS scatter of ~ 300 m s _1 per line (approximately 
one thirtieth of a pixel) . 

After discarding orders severely contaminated by sky 
emission lines, CCD fringing, or cross-disperser order 

3 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation 



overlap, as well as crowded or low-throughput orders, we 
found the useful wavelength coverage to be ^4400-6200 
A. Since the data were taken without dividing the expo- 
sures to aid cosmic ray removal (so-called "splits"), we 
developed an idl code to detect and record the locations 
of cosmic ray events in the extracted spectra. We made 
no attempt, however, to remove these or adjust the flux 
in these pixel locations but simply ignored compromised 
wavelength bins during the analysis that followed. 

4. ANALYSIS METHODS 

4.1. Classifying the cool companions 

We estimated the spectral types of the cool compan- 
ions by comparing the observed HET spectra to high- 
resoluti on templates in the ELODIE 3.1 stellar spectral 
library (jPrugniel fc Soubiranll2001h . Spectral types F0- 
M9 were considered in luminosity classes IV and V. Us- 
ing our own IDL code, we compared the best individual 
spectrum for each target to the template spectra over the 
wavelength range 5075-5200 A, which includes important 
discriminators such as the Mg 16 triplet (+MgH) and Fc 
lines. Before any comparisons were made, we degraded 
the resolution of the ELODIE templates (R=42,000) by 
convolving them with Gaussians to match the resolu- 
tion of the HET/MRS observations (R=10,000). We 
constructed a library of possible "composite-spectrum" 
models, diluting each ELODIE template spectrum by a 
flat continuum (representing the contribution from the 
hot subdwarf) using a range of dilution factors (5-95%). 
We then normalized the observed MRS spectra with 
high-order polynomials using the continuum function in 
IRAF and shifted each spectrum to its rest-frame. 

For each object-template pair we subtracted the model 
from the object and computed the sum of squared resid- 
uals (SSR); the adopted best-fitting dilution is the one 
that minimizes the SSR. We repeated this process for all 
ELODIE templates and assigned a final spectral type 
for each cool companion target corresponding to the 
global SSR minimum. To test the accuracy of our fit- 
ting technique, we also classified several twilight-sky 
spectra and a dozen RV standards with known spec- 
tral types (drawn f rom the Astronomical Almanac and 
iNidever et al.l 120021 ) . Our results matched the known 
values reasonably well, falling within two subtypes of the 
published classification 80% of the time. The largest dis- 
agreement encountered in our consistency check was 4 
subtypes. 

The second column in Table Q] lists our classifications 
for the target stars, which range from F4V to G8V. 
With 80% confidence, they are accurate to within 2 sub- 
types. They are also reason ably consist ent with previ- 
ously published values (e.g., lStark|[20 05). Although we 
report all companions as dwarfs, two of the targets, PG 
0039+049 and PG 1317+123, are matched equally well 
with subgiant templates. Classifications could not be 
made for PG 1206+165 (companion lines too weak) or 
PG 1629+081 (resolved binary, companion not on fiber). 
In these cas es, we adopt the spectral types reported by 
[Starl d200l . 

4.2. Measuring the radial velocities 

We used iraf routines to find RVs of the targets. 
Our primary strategy was to exploit the abundance of 
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Fig. 1. — Example cross-correlation of a target-sky spectrum 
pair. A single order of one of the PG 1104+243 spectra (top panel) 
is shown above the same order of a twilight-sky spectrum (middle 
panel). The He I 5876 A line from the sdB star stands out in 
the target spectrum amongst a plethora of features from the cool 
companion. The cross-correlation function (bottom panel) shows a 
reasonable match between the cool companion and solar features. 
The blue line shows the best-fitting Gaussian to the correlation 
values marked with solid points. 



spectral features from the cool companion by cross- 
correlating each spectrum against a twilight sky (i.e., so- 
lar) spectrum obtained on the same night. Before this 
was done, we created sampling region files for each spec- 
trum that defined acceptable wavelength bins to use for 
the velocity fitting. In addition to ignoring bins com- 
promised by cosmic rays and sky emission features, we 
also avoided spectral regions heavily influenced by ab- 
sorption lines from the hot subdwarf component, namely 
the Balmer lines (H/3, H7), He I lines (4472, 4713, 5015, 
5876 A) and the He II 4686 A line. Using only these 
data samples, we used the iraf continuum task to fit 
the continuum and normalize it to unity. We then cross- 
correlated the normalized spectra against the solar spec- 
tra, order-by-order, using fxcor and rvcorrect to deter- 
mine heliocentric RVs. 

Figure[T]shows an example target-sky pair and its cross 
correlation function. Most targets matched the solar 
spectrum with a correlation nearing 0.8, although sys- 
tems with companions deviating significantly from G2 V 
had lower values. Our fitting process resulted in ^14 RV 
measurements per target observation, depending on the 



exact number of echelle orders used. RV errors were com- 
puted from the standard deviation of these values about 
the mean. We find typical errors of 600-1100 m s . 

In some cases, it was also possible to estimate the sdB 
star's RV from the available He I lines and, for the hotter 
stars, He II 4686 A. We avoided using H/3 for this purpose 
since its profile spans almost an entire echelle order; also, 
its core is blended with the H/3 profile of the cool compan- 
ion. For each of these targets, after selecting sampling 
regions around the sdB He lines, we cross-correlated each 
spectrum against a self-template, chosen to be the high- 
est S/N spectrum obtained for the target. This process 
only provides relative RV measurements with respect to 
this template; converting them to absolute velocities re- 
quires determining the heliocentric velocity of the tem- 
plate. Fitting simple Gaussians or other symmetric func- 
tions to the sdB lines cannot be used for this purpose 
since the He I lines are asymmetric and would skew the 
measured centroids. Instead, we cross-correlated the ob- 
served profiles against a model sdB spectrum with T e ff 
= 30,000 K, log g = 5.6, and logN(He)/N(H) = -2. We 
then adjusted the relative measurements appropriately. 
Naturally, the RV precision achieved from this process 
is not ideal, but it was sufficient to detect and measure 
sdB accelerations in many of the systems. RVs deter- 
mined from He I lines may show a systematic offset, how- 
ever, as some of the lines are sensitive to pressure shifts 
in the high-gravity atmosphere (especially the 5876 A 
D3 line). Additionally, the gravitational redshift differ- 
ence between the cool companion and the compact hot 
subdwarf can give the appearance of a ~2 km s _1 off- 
set. These issues have be en encountered b efore in stud- 
ies of hot subdwarf s fe.g.. | Orosz et aLlll997f) and also DB 
white dwarfs (e.g., lKoesterlll987f ). where the offset can 
approach tens of km s _1 . We will discuss their effects on 
our results object-by-object, as needed. 

To investigate the zero-point accuracy of the MRS 
velocities, we cross-correlated all of the observed RV 
standards against twilight sky spectra (usually taken 
on the same night) and compared our results to the 
iNidever et ail (|2002l ) values. A histogram of the RV dif- 
ferences is described reasonably well by a Gaussian distri- 
bution with a dispersion of a k, 340 m s ; its centroid 
is within 100 m s _1 of zero. We performed a similar 
test using all of the available twilight sky spectra, cross- 
correlating them against the highest S/N twilight sky 
spectrum. In this case, we find a scatter with a — 250 m 
s _1 . Our RV standard star scatter is roughly consistent 
with the twilight sky res ult, once the addition al scatter 
of w 150 m s^ 1 from the INidever et all (|2002f ) measure- 
ments is added in quadrature. 

4.3. Fitting the orbital parameters 
We fitted each RV curve with a model of the form: 



Kad = 7 + K[cos(9 + lj) + e cos(o;)] 



(1) 



where lj is the periastron angle, K the RV semi- 
amplitude, 7 the systemic velocity, and e the eccentric- 
ity. The true anomaly (9) is defined in terms of the 
period (P), eccentricity, and time of periastron passage 
(T ) through the relation 
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tan (|) = ^tan (§ ) (2) 
and Kepler's equation 

S-esinE = M= — (t-T ). (3) 

where -E is the eccentric anomaly and M the mean 
anomaly. We used th e iDL-based rvlin packag43 
(jWright fc Howard! l2009h . which uses a Levenberg- 
Marquardt algorithm for the non-linear parts of this 
optimization problem. Orbital parameter uncertainties 
were estimated using a bootstrapping technique driven 
by rvlin (Wang et al. 2012). We determined the best 
solution with all parameters left free, and we also fitted 
the data assuming a circular orbit (e = fixed), with 
two fewer degrees of freedom. 

We also computed the floating-mean periodogram 
(e.g., ICumming et al.l [1999) for each RV data set and 
investigated \ 2 as a function of orbital period over the 
range P = 0.1-10,000 d. We consider the orbital period 
"solved" if (i) the best-fitting orbit has a \ 2 value at 
least 10 times smaller than the second-best alias and (ii) 
the uncertainty in P for the best fit is less than 5%. Once 
the first condition was met, we used a range of initial pa- 
rameter guesses in RVLIN to be certain we converged on 
the correct solution, sometimes changing the inputs only 
slightly. 

4.4. Circular or eccentric? 

Determining the eccentricity of long-period sdB bina- 
ries is an astrophysical imperative. As already mentioned 
in Sjl] if the main sequence companions are responsible for 
the formation of the sdB (and the sdB progenitor filled 
its Roche lobe), one would expect their orbits to have 
e=0 due to the circularizing nature of the RLOF and CE 
processes. If, however, the hot subdwarf formed from the 
merger of a binary originally part of a triple-star system, 
the cool companion currently observed would have had 
no part in the formation of the hot subdwarf, in which 
case the orbit could be non-circular. All of the long- 
period systems studied to date have sufficient room to ac- 
commodate an inner binary, even after al lowing for mod- 
erate "adiabatic expansion" of t he orbit (|Eggleton et al.l 
Il989t iDebes fc Sigur dsson 2002) resulting from mass loss 
in the sdB formation process. 

In order to determine whether eccentric orbit fits are 
pr eferred for our targets, w e performed an F-test similar 
to ILucv &: Sweeney! ()1971l ) using the statistic defined by 

■p (Xcirc ~ Xecc) DOF ccc . . 

~ (DOF ccc - DOF cilc ) Xe 2 cc 

where DOF ccc and DOF clrc are the degrees of freedom 
in the eccentric and circular orbit models, respectively, 
and Xecc an d Xcirc the chi-squared values of their fits to 
the data. Under the assumption that T follows an F- 
distribution, this statistic provides a quantitative means 
for deciding whether we can justify fitting the data with 
two additional degrees of freedom (e, uj) compared to our 
null hypothesis, a circular orbit solution. As T increases, 

4 http://exoplanets.org/code/ 



so does the likelihood that one can discredit the null hy- 
pothesis. We reject a circular orbit solution if the prob- 
ability of obtaining the observed J- (the p-value) falls 
below our false-rejection probability, which we conser- 
vatively set to 10 ~ 4 due to the considerable phase gaps 
in the RV measurements. If the p-value is above this 
limit, we do not rule out a circular-orbit solution using 
the current data. 

4.5. Computing the mass ratio 

For those systems for which we can measure the RVs 
of both components, we can calculate the mass ratio q 
of the binary. It is straightforward to show that the 
instantaneous sdB and main-sequence (MS) companion 
RVs at time U are related by 

V M s(U) =7(1 + 9) -qVsdBiU) (5) 

where q = K M s/K sd B = M sdB /M M s, and 7 sd B = 7ms 
has been assumed. Thus, q can be measured from the 
slope of a linear fit to a plot of RVms versus RV s dB ■ The 
resulting q is unaffected by any systematic offsets be- 
tween the component velocities, such as that caused by 
different gravitational redshifts. It is also independent 
of the inclination angle, which remains unknown for all 
systems studied (none show eclipses). We used the IDL 
fitex^ function, an o rthogonal distanc e regression rou- 
tine (procedure from iPress et al.lll992[ ). to determine q 
from the best straight-line fit to the data, taking into ac- 
count the errors from both radial velocity measurements. 

5. RESULTS 

All but two systems monitored in our program show 
some sort of orbital acceleration; most appear to have 
periods in excess of 100 days. Given such long periods 
and insufficient phase coverage, we cannot fully solve for 
the orbital parameters of the majority of targets using 
the 2005-2008 data alone. Table [2] summarizes our ini- 
tial survey results, including the maximum (absolute) RV 
change observed for each target, the range of possible or- 
bital periods, and the best-fitting orbital period taken 
from the model with the lowest x 2 in the periodogram 
(given in parentheses for the unsolved systems). Please 
note that the solutions for the unsolved systems are pre- 
liminary; they may change significantly with the addi- 
tion of new data. The four binaries previously reported 
as resolved are likely triple-star systems since their ob- 
served accelerations are incompatible with the linear sep- 
aration distances implied by their angular separations on 
the skj0. Two of these systems show RV variations of 
the main sequence star (PG 1040+234, PG 1718+519), 
while the other two show significant sdB motions (PG 
1253+284, PG 1629+081). Previous observations of PG 
1253+284 sho wing accelerat i on of the sdB but not the 
companion led iHeber et al.l (|2002| ) to a similar triple- 
star system hypothesis for this system. We are currently 
conducting follow-up observations of unsolved systems 
with HET/HRS and will report their orbital parameters 
in future publications, as they become available. 

012 began a similar RV monitoring program with Mer- 
cator/HERMES in 2009, shortly after the HET observa- 

5 http : //idlastro . gsf c . nasa . gov/ 

6 Estimates of the minimum projected linear separation dis- 
tances for the candidate triple systems range from 200 to 800 AU. 
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TABLE 2 
Preliminary Survey Results 



Target 


No. a 


RV 


Period 


Ai?.V max 




Obs. 


variable? 


[days] 


[km s" 1 ] 


Unsolved Systems 


PG 0039+049 


7 


yes" 


150-300 (220) 


11 


PG 0110+262 


9 




? 




rn 8/83 


8 








PHL 1079 


5 


yes 6 


80-500 (122) 


24 


PG 0232+095 


6 


yes 6 


<200 (14) 


29 


PG 1040+234 


9 


yes 6 


>500 (700) 


12 


PG 1206+165 


6 


yes 6 


7 


15 


PG 1253+284 


11 


yes c 


<220 (161) 


93 


PG 1449+653 


8 


yes 6 


750-1000 (859) 


15 


PG 1629+081 


8 


yes c 


2-11 (2.9) 


107 


PG 1701+359 


14 


yes 6 


630-830 (725) 


9 


PG 1718+519 


8 


yes 6 


200-330 (300) 


7 






Solved Syst 


sms 




PG 1104+243 


10 


yes 6 > c 


753.2 


16 


PG 1317+123 


11 


yes 6 - c 


1179 


13 


PG 1338+611 


11 


yes 6,c 


937.5 


16 



"number of useful HET/MRS observations from 2005-2008 
6 cool companion acceleration measured. 
c hot subdwarf acceleration measured. 



tions reported here ended. Three of their targets (PG 
1104+243, PG 1317+123, PG 1338+611) overlap with 
the HET sample. By phasing together both observation 
sets, we are able to improve upon their original system 
constraints significantly. In the case of PG 1104+243, 
we also include RVs obtained in 1995 with the Gold- 
Cam spectrograph on the 2.1m telescop e at Kitt Peak 
National Observatory (|Orosz et al.lll997l ). We discuss all 
three systems and their updated orbital parameters be- 
low. Tables [3] and |4] present the heliocentric RVs of the 
cool companion and hot subdwarf, respectively, while Ta- 
ble[S]gives the fitted system parameters described in §4.3 
along with their formal errors from the multi-parameter 
fits. 

The systemic velocity (7) is subject to a possible zero- 
point offset, so although we quote only the formal errors 
in the table, we refer the reader back to the discussion 
of our zero-point checks in ^4.21 for a more complete un- 
derstanding of their likely uncertainties. Due to the low 
quality of the hot subdwarf RVs, the orbital parameters 
shown (P, Kms, 7> e > w > ^o) were derived from the cool 
companion measurements alone. The -KsdB values re- 
ported were calculated from Kms and q, as described in 
§4.5; in all three cases, they are consistent with direct or- 
bit fits to the hot subdwarf RVs if P is taken from the cool 
companion orbital fit and the relation w s dB = wms + tt 
is imposed. Finally, we detected no significant velocity 
offsets between telescopes, so we fixed their relative zero- 
point offsets to zero during the final orbit fitting. 

5.1. PG 1104 +243 

Previous classifications for the cool compan ion of PG 
1104+243 have rang ed from F8 V (|Starll [2005') to K3.5 V 
(jAllard et all [1994). Comparisons of the HET/MRS 
spectra to the ELODIE templates show a best-fitting 
spectral type near G2 V. The dilution of the companion 
spectrum is D = L comp / L to tai ~ 0.5, where the lumi- 
nosity ratio refers to the 5075-5200 A spectral region. 



As noted by 012, the hot subdwarf in PG 1104+243 
might be described best as an 'sdOB' star since its spec- 
trum shows the He II 4686 A line in addition to strong 
H B aimer and He I featu res. Previou s RV m easurements 
have been published by lOrosz et al.l (|1997l ) , who found 
no significant variations over a three-day timespan. An 
orbital solution with P = 752 ± 14 d was recently pub- 
lished b y 012. Combining our data with those of both 
012 and lOrosz etail ((1997), we find an orbital period of 
P = 758.8 ± 0.13 d, in excellent agreement with the 012 
result. This value falls frustratingly close to two years, 
as clearly seen in the RV curve shown in Figure [2j By 
chance, current observing seasons align with quadrature 
phases; otherwise, all points would have been obtained 
near a conjunction, a less-than-ideal situation for mea- 
suring the mass ratio from spectroscopy. At the same 
time, information near conjunction is needed to improve 
constraints on the geometry of the orbit. Several more 
orbital cycles must elapse before complete phase cover- 
age can be obtained. At this time, we have no reason to 
prefer an eccentric orbit over a circular one. 

The presence of strong He lines from the sdB (5876 
& 4686 A) permits us to measure accelerations for both 
components. As expected for stars in a binary, the sdB 
and MS velocity curves are 180 degrees out of phase. The 
sdB RV curve shows a systemic velocity consistent with 
that measured from the companion (7 = — 15.5±0.03 km 
s -1 ). From the relationship between RV s dB and RVms, 
as seen in Figure^ we compute a mass ratio of q = 0.68 
± 0.08. Table [5] summarizes the system parameters de- 
rived for PG 1104+243. If we assume a companion mass 
of 1.0 Mq and adopt a 10% error to include uncertain- 
ties in the spectral classification, we find M s dB = 0.68 ± 
0.11 Mq, considerably larger than the 'canonical' value 
of ~0.48 M Q , which corresponds to helium core ignition 
at the tip of the RGB. If we instead assume exactly the 
canonical mass for the hot subdwarf, we would compute 
a mass of 0.7 ± 0.1 Mq for the companion, which implies 
its spectral type might be significantly later than our re- 
sult from the composite spectrum analysis, or that the 
companion is evolved. 

5.2. PG 1317+123 (Feige 80) 

Spectral classifications for the hot component in PG 
1317+123 have varied over years. It is probably an sdO 
star as it displays a strong He II 4686 A line and a rela- 
tively weak He I 4472 A line. lUlla fc TheiU (fl99l clas- 
sify the cool companion as a G8V, while 012 report it 
to be a 'G' star. We find a spectral type near Gl V and 
dilution D w 0.3 from our comparisons with ELODIE 
models. 012 reported an orbital period of P — 912 ± 
52 d. After phasing the HET data with theirs, we find 
a significantly longer value, P = 1179 ± 12 d. Figure 
[5] shows the resulting RV curve, while Table [5] summa- 
rizes all of the orbital parameters derived. Our results 
for the systemic velocity and MS RV amplitudes also dis- 
agree. As pointed out by 012, the Mercator/HERMES 
measurements show much larger scatter than one would 
expect given the precision of their data. They briefly 
suggested 7 Doradus-type pulsations or star spots on 
the cool companion as potential culprits. However, G- 
type classifications for the cool companion are inconsis- 
tent with those of 7 Dor stars, which are typically early F 
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TABLE 3 

Cool Companion Velocities for Solved Systems 



HJD 


RVms 


Facility" 


HJD 


RVms 


Facility" 


HJD 


RVms 


Facility" 


-2450000 


[km s _1 ] 




-2450000 


[km s _1 ] 




-2450000 


[km s _1 ] 





PG 1104+243 PG 1317+123 PG 1338+611 



-215 


23370 


- 1 6 


SO 


_|- 


4 


70 


K 


3379 


93771 


45 


00 




1 


09 


H 


3379 


ou,k) 


24 


89 


_|- 


Q 


90 


II 


-214. 


16710 


-19 


.30 


-|- 


4 


.50 




3392 


92593 


44. 


.28 


± 


x. 


00 


H 


3400 


89553 


25 


.(>2 


-|- 


x. 


21 


H 


-213. 


24680 


-17 


.30 


-|- 


4 


.70 


K 


3421 


.84076 


44 


.04 


± 


0. 


68 


H 


3415 


.89117 


27. 


.84 


-|- 


0. 


77 


H 


-212 


32400 


-17 


60 


-|- 


4 


.40 


K 


3446 


76922 


44 


01 


± 





63 


H 


3429 


87462 


98 


35 


-|- 





57 


H 


3379 


84221 


-20 


53 


-|- 


o 


.74 


H 


3476 


68120 


42 


92 


± 





99 


H 


3476 


70485 


31 


35 


-|- 





63 


H 


3302 


80703 


-20 


20 


-|- 


o 


.72 


H 


.'5488 


65464 


42 


87 







91 


H 


34S 8 


6698(1 

• U U ir O 17 


39 


69 


-|- 





96 


H 


3415 


95642 


-20 


06 


-|- 


o 


56 


H 


3500 


62543 


42 


13 


± 





70 


H 


3502 


64537 


39 


.77 


-|- 





64 


H 


34 9 9 


73480 


-19 


63 


-|- 


o 


.41 


H 


3724 


01307 


36 


39 


± 





51 


H 


3731 


99994 


40 


93 


-|- 





78 


H 


3447 


65376 


-18 


90 


-|- 


o 


54 


H 


3771 


87931 


35 


.42 


± 





78 


H 


3732 


97699 


40 


65 


-|- 





73 


H 


3462 


69^99 


-19 


61 


-|- 


o 


.40 


H 


4453 


00102 


48 


53 


± 


I 


23 


H 


4479 


96424 


35 


36 


-|- 





78 


H 


3476 


.78809 


-19 


.43 


-|- 





.51 


H 


4551 


.75800 


45 


.21 




0. 


90 


H 


4540 


.80436 


37 


03 


-|- 


0. 


83 


H 


3498 


.74435 


-19 


.12 


-|- 





.46 


H 


5030 


42431 


34 


.08 




0. 


.09 


M 


5028 


.44884 


27. 


83 


-|- 


0. 


02 


M 


3757 


8966^ 


-11 


03 


-|- 


o 


.40 


H 


5030 


42431 


34 


08 


± 





09 


M 


5028 


47026 


27 


.77 


-|- 





09 


M 


4479 


82452 


-10 


.97 


-|- 


o 


.78 


H 


5030 


45097 


39 


69 







17 


M 


5028 


49174 


27 


.77 


-|- 





03 


M 


5204 


70741 


-12 


24 


-|- 


o 


09 


M 


5030 


45097 


39 


69 


± 





17 


M 


5341 


57223 


30 


81 


-|- 





03 


M 


5204 


72527 


-12 


.37 


-|- 





.04 


M 


5294 


.54985 


39. 


10 


-j- 


0. 


10 


M 


5351 


.41840 


3 1 


.07 


-|- 


0. 


.05 


M 


5217 


.67376 


-12 


.29 


-|- 





10 


M 


5294 


.54985 


39. 


10 


-j- 


0. 


10 


M 


5371 


.41748 


32 


84 


-|- 


0. 


(12 


M 


5217 


69174 


-12 


.27 


-|- 


o 


10 


M 


5353 


40135 


39 


53 


± 





1 1 


M 


5566 


78096 


40 


00 


-|- 





04 


M 


5234 


60506 


-11 


.42 


-|- 





.04 


M 


5353 


.40135 


39. 


.53 


-j- 


0. 


1 1 


M 


5616 


64021 


40 


.29 


-|- 


0. 


02 


M 


5234 


.61608 


-12 


.09 


-|- 





.11 


M 


5373 


.41275 


41. 


12 


-j- 


0. 


08 


M 


5616 


64021 


40 


.29 


-|- 


0. 


02 


M 


5234 


62710 


-12 


02 


-|- 


() 


12 


M 


5373 


41275 


41 


12 


-j- 





OS 


M 


5622 


70741 


39 


98 


-|- 





04 


M 


5264 


.54423 


-11 


.15 


-|- 





.07 


M 


5566 


.74478 


46. 


66 


-j- 


0. 


1 3 


M 


5622 


73169 


39. 


96 


-|- 


0. 


03 


M 


5340 


42821 


-11 


.93 


-|- 





.07 


M 


5566 


.74478 


46. 


66 


± 


0. 


1 3 


M 


5656 


.60318 


39. 


.64 


-|- 


0. 


04 


M 


5351 


37760 


-11 


91 


-|- 


() 


08 


M 


5616 


61438 


47 


.49 


-j- 





06 


M 


5658 


58702 


39 


98 


-|- 





OS 


M 


5553 


70059 


-18 


20 


-|- 


() 


1 1 


M 


5616 


61438 


47. 


.49 


-j- 





06 


M 


5660 


61069 


39 


86 


-|- 





04 


M 


5569 


69768 


-18 


.78 


-|- 


n 


09 


M 


5616 


61438 


47 


.49 


-j- 





06 


M 


5663 


54111 


39 


90 


-|- 





03 


M 


oo / y 


. OOOO..J 


1 Q 

-_i_y 


99 


_i_ 


n 


1 2 


M 


OOlU 


.Olloo 


ft 1 . 


. ^y 




n 
u. 


(\a 
uo 


M 


0000 


.OO^iOO 


'i G 

■ >)) 


. 00 




n 
u. 


(19 


M 


5579 


.58348 


-18 


.93 


± 





.09 


M 


5640 


.61157 


45. 


.88 


± 


0. 


06 


M 


















5589 


.75559 


-19 


.06 


± 


0. 


.10 


M 


5640 


.61157 


45. 


.88 


± 


0. 


06 


M 


















5611 


.64081 


-19 


.73 


± 





.07 


M 


5653 


.60091 


44 


.58 


± 


0. 


11 


M 


















5622 


.61268 


-19 


.95 


± 





.08 


M 


5653 


.60091 


44 


.58 


± 


0. 


11 


M 


















5639 


.58778 


-19 


.57 


± 


0. 


.07 


M 


5661 


.63454 


45. 


.09 


± 


0. 


09 


M 


















5650 


.50737 


-20 


.16 


± 





.07 


M 


5661 


.63454 


45. 


.09 


± 


0. 


09 


M 


















5652 


.55550 


-19 


.99 


± 


0. 


05 


M 


5722 


.45431 


44 


.93 


± 


0. 


24 


M 


















5655 


.60537 


-19 


.87 


± 


0. 


.07 


M 


5722 


.45431 


44 


.93 


± 


0. 


24 


M 


















5658 


.47913 


-20 


.10 


± 





.06 


M 


































5664 


.51354 


-20 


.17 


± 





.07 


M 


































5666 


.43280 


-20 


.41 


± 





.05 


M 


































5666 


.45153 


-20 


.32 


± 





.06 


M 


































5672 


.53688 


-20 


.30 


± 





.11 


M 


































5685 


.44968 


-20 


.02 


± 





.06 


M 


































5718 


.44565 


-19 


.50 


± 





.06 


M 



































a K = KPNO 2.1m/GoldCam spectrograph fOrosz et al. 1997 V H=HET/MRS, M=Mercator /HERMES (012). 



or late A dwarfs. The other scenario, a bright spot on the 
companion's surface, seems more compatible with the ob- 
servations. Assuming a solar-type dwa rf and a star spot 
with a filling factor of a few percent ()Saar fc Donahue] 
Il997f ). the amplitude of the scatter (1-2 km s _1 ) would 
imply a rotational period around one week. As with 
PG 1104+243, we currently have no reason to prefer an 
eccentric-orbit solution. Even with additional data, con- 
straining the eccentricity will be especially difficult until 
the rapid velocity variations are better characterized. 

He I 5876 A and He II 4686 A features provided RV 
estimates for the hot subdwarf that are again 180 degrees 
out of phase with the cool companion. Figure [3] shows 
a strong linear relationship between KV sc io and RVms, 
whose slope gives q = 0.42 ± 0.05. A circular orbit fit- 
ted to the sdO RVs with the period and phase fixed to 
the values in Table [5] gives a consistent systemic velocity 



and a semi-amplitude in agreement with the expecta- 
tion from combining Kms and q. Assuming a companion 
mass from with our spectral classification, Mms = 1-03 
± 0.11 Mq, gives M sd0 = 0.43 ± 0.07 M Q . 

5.3. PG 1338+061 (Feige 87) 

Our analysis of PG 1338+061's spectrum reveals a G2- 
G7V companion (best fit: G4V with D w 2) much 
earlier than previou s classifications fKOV. [StarH [20051 
K4.5V. lAliard et al.l[l99l . Two He I lines (5876 & 4472 
A) were used to measure the sdB RVs, which show the 
sdB and companion accelerating in anti-phase with one 
another. The last panel in Figure[2]presents the RV curve 
for PG 1338+061, combining the HET and Mercator ob- 
servations. Unlike PG 1104+243 and PG 1317+123, we 
find significant evidence for an orbit with non-zero eccen- 
tricity. Figure [4] shows our best-fitting circular and ec- 
centric fits to the cool companion's RV curve for compar- 
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TABLE 4 

Hot Subdwarf Velocities for Solved Systems from HET/MRS 



HJD RV MS HJD RV MS HJD RV MS 

-2450000 [km s" 1 ] -2450000 [km s" 1 ] -2450000 [km s" 1 ] 

PG 1104+243 PG 1317+123 PG 1338+611 



3379.84221 


-7.8 ± 1.9 


3379.93771 


28.8 


± 


1.3 


3379.95036 


41.8 


± 


3.6 


3392.80703 


-8.0 ± 1.6 


3392.92593 


31.9 


± 


2.6 


3400.89553 


43.3 


± 


3.3 


3415.95642 


-9.6 ± 2.8 


3421.84076 


31.4 


± 


1.0 


3415.89117 


39.1 


± 


2.9 


3429.73480 


-11.1 ± 2.3 


3446.76922 


32.2 


± 


2.2 


3429.87462 


38.3 


± 


3.6 


3447.65376 


-11.3 ± 1.1 


3476.68120 


31.8 


± 


1.0 


3476.70485 


32.5 


± 


2.4 


3462.62302 


-11.4 ± 2.0 


3488.65464 


34.6 


± 


3.0 


3488.66280 


31.4 


± 


3.5 


3476.78809 


-12.0 ± 2.3 


3500.62543 


37.6 


± 


2.9 


3502.64537 


33.1 


± 


4.4 


3498.74435 


-11.8 ± 2.4 


3724.01307 


50.1 


± 


1.8 


3731.99294 


18.4 


± 


3.6 


3757.88245 


-22.5 ± 1.5 


3771.87931 


51.3 


± 


5.4 


3732.97699 


16.0 


± 


2.1 


4479.87074 


-23.6 ± 2.3 


4453.00102 


24.2 


± 


8.3 


4479.96424 


28.3 


± 


3.9 






4551.75800 


31.2 


± 


3.0 


4540.80436 


20.1 


± 


6.2 



TABLE 5 

Orbital Parameters for the Solved Systems 



Target 


P 


Kms 




7 


e 






T C 


Name 


[days] 


[km s _1 ] 


[km s" 1 ] 


[km s _1 ] 




[dcg] 


[HJD 


-2450000] 


PG 1104+243 


753.2 ± 0.8 


4.43 ± 0.06 


6.5 ± 0.8 


-15.68 ± 0.05 


0.0" 




4534.5 ± 3.3 




PG 1317+123 


1179 ± 12 


6.2 ± 0.2 


15.5 ± 1.7 


+40.3 ± 0.2 


Q.0 b 




4453 ±15 




PG 1338+611 


937.5 ± 1.1 


8.8 ± 0.15 


15.2 ± 1.8 


+32.58 ± 0.07 


0.15 ± 0.02 


209 ± 6 


4685.5 ± 3.0 





"computed from Kms an d <?■ 

'fixed; no reason to prefer eccentric fit over circular solution at this time. 

c time of velocity maximum (maximum redshift) in the cool companion RV curve (e = case) or time of periastron passage (e ^ 0). 



ison. Computation of the F-test statistic gives J-~=12.6 
with a corresponding p-value of 10 -7 , well below our cri- 
terion for claiming non-zero eccentricity. We therefore 
reject the circular-orbit hypothesis and report e = 0.15 
± 0.02 for this system. The best-fitting period in this 
case is P = 937 ± 10 d, in slight disagreement with 012, 
who assumed a circular-orbit solution and found P = 
915 ± 16 d. In t j6.2l we discuss the implications of non- 
zero eccentricity on formation scenarios of long-period 
sdB systems. 

The relationship between RV s< \b and RVms (Figured]) 
is fitted well with q = 0.58 ± 0.07. Combining this value 
with Kms predicts a velocity semi-amplitude consistent 
with the best-fit to the subdwarf RV curve (with period 
and phase fixed). The difference in the systemic veloc- 
ities is 7 s dB — 7ms = — 2.1± 1.0 km s . Assuming a 
mass of Mms = 0.97 ± 0.1 Mq for the companion, we 
calculate Af s dB = 0.56 ± 0.09 M© for the hot subdwarf, 
consistent with the canonical value. 



6. THE PERIOD DISTRIBUTION OF 
HOT SUBDWARF BINARIES REVISITED 

Recent orbital solutions for long-period binaries have 
drastically changed the empirical period distribution of 
hot subdwarf binaries, compared with the situation a 
few years ago. We present in Figure [5] an updated 
period histogram for all systems with published or- 
bital parameters. Most of t h e shor t-period systems 
can be found in iGeier et all ft20lli) ( see Table A.l 
and references therein) . iCopperwheat "etaLl (Hoiil), and 



iMorales-Rueda et al.l ((2003). For the long-period b ina- 
ries, which include PG 1018-047 (|Deca et al.ll2012T) . re- 
sults from 012, and all targets in our survey with rea- 
sonably well-constrained periods (see Table [2]), we in- 
clude estimates of the uncertainties in the orbital peri- 
ods. Given the difficulties associated with monitoring 
the long-period systems, several strong observational bi- 
ases are undoubtedly at play in the diagram; premature 
conclusions should not be drawn. 



6.1. The problem posed by long orbital periods 

One clear trend appears to emerge from the empirical 
distribution of orbital periods: closer binaries tend to 
have late-type MS or WD companions while the long- 
period systems host F/G/K dwarfs. This appare nt di- 
chotomy has been noted before fe.g.. IGreen 1 120081 ) . Al- 
though the number of long-period binaries known is still 
relatively small, their addition allows us to begin test- 
ing the predictions of BPS models in both the short- and 
long-p eriod regim e s for the firs t time. 

The lHan et all (|2002l . |2003[ ) models reproduce the 
short-period end of the distribution quite well but pre- 
dict that the sdB+MS (G/K) binaries should have peri- 
ods shorter than ^100 d, the outcome of CE evolution 
using the "a-formalism" , in which orbital energy (and/or 
thermal energy stored in the envelope) is tapped to ex- 
pel the envelope. Current observations suggest that the 
energy-based a-formalism for CE evolution in its sim- 



9 



-12 



> 



1 1 1 

l'\ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i — ' 

l'\ l\ i'\ f\ I 


1 1 1 1 1 1 1 1 1 1 1 1 1 r~ 

PG 1104+243 
t * ^ '' ' „ „ ' 

f A.' i A' [ f\i - A/ 


■ V 
A j 


A .' 1 A .' 1 A ' 1 A ' 1 A /' 

/ V '/ V / V / V y V 
\ \ k \ k \ k \ k \ 


/ '\ / A / '\ / h 


" ' VJ 


* / vy ^ / \j 1 / v/ * / \j ^ i\ 
\i \j \j \j \j * 


m \ i \ \ ■ \ / 1 ■ I J 1 1 




— i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i 


1 1 1 1 l-M 1 1 1 1 1 1 1 V~ 






i « • c. * 


B 


.- -i 

i .... i .... i .... i 





1000 2000 3000 4000 5000 6000 




^ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1- 



3500 



4000 4500 5000 

Time [HJD - 2450000] 



5500 



Fig. 2. — Heliocentric RV curves for the three targets overlapping with the study of 12: PG 1104+243 ( top), PG 1317+123 (middle), 
and PG 1338+611 (bottom). Cool companion measurements from HET/MRS, 012, and Orosz ct al. ( 1997]) are shown with filled circles, 
squares, and triangles, respectively. Open circles represent sdB velocities measured from He lines in the HET/MRS spectra. The best- 
fitting orbital solutions for the sdB and cool companion are denoted with dot-dashed (orange) and solid (blue) lines, respectively. Residuals 
from this fit are shown in the lower portion of each panel for the cool companion measurements. Error bars on the Mercator/HERMES 
points are much smaller than the symbol sizes; we do not plot them here. 



plest forms does not ac count for the long-peri od sdB 
binaries. We note that lHan et all (|200a 12003ft delib- 
erately excluded the GK-effect binaries when choosing 
their favored BPS models, focusing on the short-period 
M dwarf and WD systems (the only systems with mea- 
sured periods available at the time). Thus, it is not 
surprising their model predictions are inconsistent with 
the c urrent period distribution of sdB+F/G/KV bina- 
ries. I Clausen et al.l (|2012[ ) encountered a similar prob- 
lem when modeling the sdB population using the a— 
formalism. This failing is not limited to the sd B binaries 
and h as been discussed in other contexts by iWebbinkl 
(2008), who points out the need in some cases for (1) 



additional energy sources to augment or substitute for 
some of the energy extracted from the orbit in expelling 
the envelope, or (2) a reduction in the binding energy of 
the envelope (possibly through more extensive mass loss 
prior to the CE phase), in order that orbital shrinkage 
not be so drastic. 

An altern ative scenario proposed by iNelemans et al.l 
(|2000l 120011 ) appears, at first glance, to be more suc- 
cessful. In this model, the long-period systems are the 
product of CE evolution, with a prescribed fractional an- 
gular momentum loss that is proportional to the frac- 
tional mass lost when the envelope is ejected. The cur- 
rent distribution of long-period systems overlaps surpris- 
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Fig. 3.— HET/MRS orbital velocities of the cool companion and sdB plotted against one another for PG 1104+243, PG 1317+123, and 
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Fig. 4. — Heliocentric radial velocities for PG 1338+611 with the 
best-fitting circular (blue dotted line) and eccentric (solid black 
line) orbits, and their residuals (bottom panels). Velocities from 
HET/MRS and 012 are shown with circles and squares, respec- 
tively. The data are not consistent with a circular solution; the 
best-fitting eccentric model has e = 0.15 ± 0.02. 



ingly well with the densest region of the predicted pop- 
ulation whe n the angular mo mentum parameter 7 = 1.5 
(Figure 2 of INelem ans 2010). Current observations hint 
at a dearth of binaries in the range P ~ 30-110 d, and 
in this context we note that a period histogram drawn 
'by eye' from the sdB+MS po pulation (with co mpanion 
masses ~ 1-&T©) predicted by INelemans! (|2010f ) shows a 
local minimum in the number of binaries with periods 
from « 30 to 90 d. It would be well not to read too 
much into this apparent agreement, since the model is 
advertised as tentative, is not well-matched to the PG 
sample in terms of its flux limit, and the transformation 
of sdB bolometric luminosities to V magnitudes is not 
described. In addition, the empirical diagram should be 
more thoroughly populated and the observational biases 
accounted for, before the significance of the gap is as- 
sessed. 




Log P [days] 

Fig. 5. — Histogram of the orbital periods of 140 hot subdwarf 
binaries with known orbital parameters. The black vertical lines 
at the top indicate the individual periods of all systems. For the 
long-period binaries with poor fits, they represent the best-fitting 
orbit in the periodogram; two diagonal lines emanating from the 
tops of these symbols mark the range of possible orbital periods. 
In the rows of vertical lines below the first, we re— mark systems 
with F-K dwarf (second row, orange marks), M dwarf (third row, 
red marks), and white dwarf (fourth row, blue marks) companions. 
Those with uncertain companion classifications or none at all are 
not shown. The long-period systems exclusively have companions 
of spectral type ~K4 or earlier. Although the statistics are poor, 
observations currently show a gap centered near P = 70 d. 



6.2. Eccentricity as a clue to the long periods 

It is interesting to consider the question of eccentric- 
ity of the observed wide orbits. PG 1338+061 (th is pa- 
per) and possibly PG 1018-047 (]Deca et al.ll2012[) have 
significant eccentricities. We also note that eccentric or- 
bits have not been excluded for PG 1104+243 and PG 
1317+123 based on the present data. If wide sdB binaries 
typically prove to have eccentric orbits, this would sug- 
gest one of three possibilities. Perhaps (1) these systems 
are or were hierarchical triple systems in which the sdB 
star was made in the inner bin ary, possibly by merger 
(see, e.g.. lClausen fc Waddl201ll who describe a "hydro- 
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gen merger channel" that leads to sdB's with canonical 
masses). We note that in PG 1338+061, the minimum 
periastron distance of 2 AU in the current orbit gives 
ample room to have once accommodated an inner sys- 
tem. If however the systems have always been binaries, 
then either (2) circularization was not achieved because 
the sdB progenitor never filled its Roche lobe or (3) ec- 
centricity was pumped into the system via one of several 
possible mechanisms. 

In the first case, the current long period eccentric or- 
bits are not directly comparable to the predictions of 
BPS models since the observed cool companions did not 
participate directly in the formation of the hot subwarfs. 
They nevertheless provide some dynamical constraints 
on "binary" formation scenarios, since they limit the 
size of the inner orbit in the hierarchical triple. Studies 
of the sdB rotational velocities in these systems might 
shed light on whether a merg er included two He-core 
WDs (fast rotation expected: iGourgouliatos fc Jeffervl 
2006; although if there is a giant stage between merger 
and the sdB phase t he rotation may be slowed, see 
iZhang fe Jeffervll2012ft or a He WD and a MS star (slow 
rotation, owing to angular momentum s hedding near the 
tip of the RGB; iClausen fc W^IMTTh . 

In the case these systems have always been binaries, we 
suppose that a st rong, tidally-enhanced w ind along the 
lines suggested bv lTout fc Eggletonl (1988) can occasion- 
ally result in the complete loss of the envelope, just as 
the core of the sdB progenitor rea ches the critical mass 
(or near it, see iD'Cruz et aT1ll996ft to ignite He burning, 
without ever filling its Roche lobe. Some expansion of 
the orbit would occur owing to the wind mass loss, and 
the orbital eccentricity could in large part be preserved. 
If this timing coincidence does not occur, the extra mass 
loss can at least reduce the envelope binding energy and, 
after a CE episode, leave a (circular) orbit at longer pe- 
riod than would otherwise be attained, or even reduce 
the mass ratio to the point where RLOF mass transfer 
pr oceeds stably and a CE phase is avoided altogether, as 
in lTout fc Eggletonl (fll)88ft . 

Finally, we consider briefly the case that these systems 
are in fact binaries in which the sdB progenitor once filled 
its Roche lobe. Several mechanisms exist that could work 
to increase the orbital eccentricity, even as mass is being 
lost or transferred, thereby preventing the binaries from 
circularizing. Most notabl y, these include pum ping from 
a circumbinary disk (e.g. iDermine et al.l 120121 and ref- 
erenc es therein) and asymmetric /ph a se-dependen t mass 
loss (jBonacic Marinovic et all 120081: ISokerl 12000ft . both 
of which have been invoked to explain the large ec- 
centricities observed in some post-AGB systems hav- 
ing undergone strong tidal interactions. Long-period 
sdB binaries are post-RGB systems, and as such, their 
sdB progenitors had significantly different luminosities, 
radii, envelope structures, and mass-loss rates than AGB 
stars. Eccentricity-pumping mechanisms operating on 
RGB and AGB systems will neither act on the same 
timescales, nor will they have the same ability to com- 
pete with the effects of tides. We regard these scenarios, 
which have primarily been linked to post-AGB binaries, 
as uncertain in the case of hot subdwarf binaries; further 
discussion of them is beyond the scope of this paper. For 
completeness, we also recognize th e action of a distant 
third body (the Kozai mechanism, |KozaiHl962D . which 



we regard as the 'ultimate epicycle', as another possible 
avenue for increasing the eccentricity. 

7. CONCLUDING REMARKS 

We have described an ongoing program to determine 
the orbital parameters of hot subdwarfs in binaries with 
F-K type main sequence companions. Preliminary re- 
sults show most of our targets have orbital periods on 
the order of years (rather than days), in agr eement with 
the re cent study by 012 and the findings of iGreen et al.l 
(2001) from more than a decade ago. Precise orbital so- 
lutions cannot be provided for most of our systems at 
this time, owing to insufficient phase coverage. Nonethe- 
less, the preliminary findings imply that the assumptions 
inherent in some subdwarf formation scenarios should be 
revised, which probably requires both re-tuning the pa- 
rameters used in BPS models and accounting for at least 
some of the long-period systems as former hierarchical 
triples. 

For the three systems also observed by 012, we calcu- 
lated precise orbital parameters from nearly seven years 
of RV measurements. Mass ratios were derived from the 
orbital RVs and are consistent with our spectral classifi- 
cations for the cool companion and canonical-mass hot 
subdwarfs for two of the three systems. In the case of 
the third, PG 1104+243, the derived subdwarf mass ex- 
ceeds the range classically expected for sdBs that evolved 
through the stable RLOF channel, implying the subdwarf 
either formed from a WD merger (resulting in a larger- 
than-canonical mass) , or the companion has a mass dif- 
ferent from what we infer from its spectral type. 

A long-term goal for observers is to determine the joint 
distribution of orbital periods and eccentricities in long- 
period sdB binary systems. At this time, we cannot rule 
out eccentric orbit solutions for PG 1104+243 or PG 
1317+132. Determining the eccentricity of PG 1317+123 
will prove especially difficult due to the presence of rapid 
RV variations, which might originate from star spots on 
the cool companion. We find a statistically significant 
eccentricity of e = 0.15 ± 0.02 for PG 1338+061 and re- 
ject the circular-orbit hypothesis in this case. The long- 
period binary PG 1018 —047 also shows ev idence in favor 
of a non-circular orbit (' Deca et al.l l2012). although this 
result needs confirmation. Such departures from e = 0, 
taken at face value, challenge the prevailing notion that 
strong binary interactions, involving the filling of the sdB 
progenitor's Roche lobe, are always involved in the for- 
mation of hot subdwarf stars. If each of these binaries 
was (or is) a hierarchical triple system, then the cool com- 
panion that we observe probably did not participate di- 
rectly in the formation of the hot subdwarf, in which case 
its orbit would not have been circularized. The hot sub- 
dwarf in this picture could have formed fro m the merger 
of two He-core WDs (lHan et al.ll2002l 12003ft o r a He-core 
WD and a MS star |ciausen fc Wadd 12011ft . Measure- 
ments of the rotational velocities of sdBs in long-period 
systems might shed light on their formation by identify- 
ing the type of merger. On the other hand, if the sys- 
tem h as always been a binary , a strong, tidally-enhanced 
wind (jTout fc Eggletonll988ft could have assisted the sdB 
progenitor in ejecting its envelope on the RGB before 
it expanded enough to fill its Roche lobe, thereby pre- 
serving some level of eccentricity in the binary. Even 
in the case the Roche lobe was filled previously, non- 
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steady mass loss, the action of a distant third companion, 
or interactions with a circumbinary disc could have po- 
tentially added eccentricity to the orbit, although these 
scenarios seem less likely. 

In light of the updated period histogram for sdB bina- 
ries (Fig. [5]), we find it imperative that observers con- 
tinue to study both short- and long-period systems with 
all possible types of companions, to arrive at a clear un- 
derstanding of the formation channels for hot subdwarf 
stars. Thus it is important to give attention to issues of 
completeness and bias in our knowledge of binary orbits 
(or otherwise) for sdB stars. An important first step is to 
establish a more complete catalogu e or finding lis t , as fo r 
example in the effort described by iGirven et al.l |2012), 
who search for photometrically composite objects and 
use decomposition of the spectral energy distribution to 
identify binaries containing h ot subdwarfs. We no te that 
one of the selection criteria of lGirven et al.1 ()2012[ ) is that 
the binary is unresolved; this by itself does not indicate 
whether the binary is close enough to be tidally interact- 
ing, so follow-up studies are needed. Eventually the goal 
is to derive full orbital elements for an unbiased sample of 
any interacting binaries that are found in these surveys. 



Meanwhile, work continues to characterize the orbits of 
sdB binaries that are already known. We are currently 
conducting a follow-up program with HET/HRS to de- 
termine the orbital periods, velocities, and eccentricities 
of the unsolved targets in our survey. 
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E. Eberly. This work is also based on observations made 
with the Mercator Telescope, operated on the island of 
La Palma by the Flemish Community, at the Spanish Ob- 
servatorio del Roque de los Muchachos of the Instituto 
de Astroffsica de Canarias. Finally, we thank an anony- 
mous referee for comments and suggestions that helped 
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